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Abstract 

A grain boundary which separates domains with different orientation along 
the c-axis is analyzed. The coupling of two parallel superconducting planes 
whose order parameters are rotated leads to interesting properties which are 
expected to be common to other kinds of grain boundary. The density of 
states is enhanced at low energies. The number of localized states depends on 
the degree of misfit across the boundary. A continuum of zero energy states 
appears at certain values of the angle of rotation. 
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High temperature superconductors have been deeply studied since their discovery in 1986 
These new materials are based on planar copper oxide structures that are insulating 
and present a phase transition to a superconducting state when they are doped. The micro- 
scopic models proposed to explain the appearance of superconductivity in these materials 
point to a non conventional symmetry of the order parameter. Now, it seems that finally 
the d-wave model is fully accepted as the correct one, at least for hole-doped super- 



conductors. The determination of the pairing symmetry has been possible thanks to the 
experiments done using different techniques. Those experiments done in samples with grain 
boundaries have been especially important. It is well known that lattice distortions, such as 
grain boundaries, provide an intrinsic source of frustration in anisotropic superconductors 
because the order parameter follows the orientation of the lattice axes. These systems have 



been widely studied in the literature both experimentally |^-|T2|,|T5|-|Tj| and p8|-p6[. First 
experimental results were controversial. Experiments reported by Sun and collaborators 
|T5fl showed the observation of a small but finite Josephson coupling between a conventional 
s-wave superconductor and a heavily twinned film of YBCO. This was interpreted as an ev- 
idence against d-wave symmetry. Josephson coupling across interfaces between two YBCO 
regions with crystalline axis misoriented by 45° led Chaudhari and Lin ]T7|] to the same 
conclusion. On the other hand, the detection of a half quantum of flux on tricrystal ring 
experiment P,|T0|| gave a clear evidence of the existence of nodes and lobes in the gap, sup- 



porting the d-wave model. Recently it has been reported an experiment []Tl| which, based 
on symmetry considerations, demostrate a pure d-wave order parameter in the tetragonal 
superconductor T12201. 

However, the orthorhombic structure of YBCO allows for the existence of d + s pairing 
symmetry. Direct evidence of this behavior has been found on the magnetic interference 
pattern of a sample with a single twin boundary |12|| , at which the direction of a and b 



axes reverse. As the junction is formed in the superconducting plane the whole system lacks 
periodicity. This fact introduces complexity in the study of the excitation spectrum. Latter 
has been studied |||] within the framework of the Bogoliubov-de Gennes (BdG) equations 



and from ad hoc microscopic models. The excitation spectrum presents a zero-energy peak 
in the local density of states near the twinning plane. Here we show that a finite density of 
states at E — can also be obtained in another kind of grain boundaries. 

In the present article, a new type of grain boundary is analyzed, see Fig. 1. We consider 
two parallel superconducting planes with d-wave pairing symmetry. One of the planes has 
its order parameter rotated with respect to the other by an angle 0o- Both planes are 
coupled by a hopping term t, which allows the exchange of electrons between the two planes. 
Interplane tunneling of Cooper pairs J27[ could also be included in the model system but, 
for the sake of simplicity, we have not taken it into account. Notice an important difference 
between our system, where the grain boundary is formed between two planes, and grain 
boundaries usually considered in the literature with only one plane involved. Due to the 
symmetry of our system, the spin and momentum are conserved in the hopping process, 
contrary to what happens in other kind of grain boundaries, where translational invariance 
is violated. This fact opens us the possibility to obtain the whole spectrum and to study 



the characteristic features of grain boundaries avoiding quasiclassical approximations |28 
Despite of its simplicity, we expect the spectrum of our system to reproduce the same 
behavior than in plane grain boundaries. Besides, this kind of structure could not only be 
found in YBCO but also in other compounds. In particular, the case <p = | is equivalent to 



7r-junctions p0|j2^ , |2^| . It has also been suggested the possibility that this type of structures 
is formed in granular materials as well as during the growth process of the superconducting 
planes. A more realistic model should consider the coupling of two semiinfinite media but 
we expect that the main properties are captured by the present model. 

We start from the mean field Hamiltonian 

ka ka 

+ ^2 ^ka c kaa C ^-k-aa + ^ka C -k-aaCkaa 
k k 



+ 2 tc La C kaf3 + Yl tC Lf3 C kaa (1) 



ka ka 



where greek indices label the planes, c\. aa creates an electron with momentum k and spin 



a in the plane a, and = — jj, with jj, being the chemical potential. This dispersion 
relation has been chosen because of its simplicity but the generalization to one with a 
different dependence on k is straightforward. The angle of rotation O comes into the 
problem through the relative values of the order parameters: 

A ka = A cos 20 

Aj^ = A cos 2(0 + 4> ) (2) 

To calculate the eigenvalues of this Hamiltonian the BdG equations are solved. After a 
straightforward calculation we obtain: 

El = i 2 k + t 2 + ^ [cos 2 20 + cos 2 2(0 + O )] 
A 4 

± { — [cos 2 20 - cos 2 2(0 + O )] 2 + 4^t 2 

+ t 2 A 2 [cos20-cos2(0 + o )] 2 }^ (3) 

This spectrum presents states at the Fermi level. These zero energy states, or nodes, are 
given by the equation: 

A 2 

H = t 2 - — [cos 2 20 + cos 2 2(0 + O )] 



± A | cos 20 + cos 2(0 + O ) 



A 2 

x J — [cos 20 - cos 2(0 + O )] 2 - t 2 (4) 

From this equation it is clear that nodes can only exist if one of the following conditions is 
satisfied: 

2t 

| cos 20 — cos 2(0 + O ) I > I — I (5) 



or 



A 



cos20 + cos2(0 + o ) = (6) 



Two parameters are relevant: the hopping t between the planes and the angle O . For small 
values of t the most relevant condition is (|5|) except for O = tt/2. On the other hand, for 
sufficiently large values of t, (|5|) is never fulfilled and the nodes are given by 

4 2 = t 2 - A 2 cos 2 20 (7) 

with given by (j6|). The solution of eq.(|]) in the first Brillouin zone is plotted in Figs. 2a 
and 2b for 0o — tt/6 where the influence of the hopping can be seen. The nodes appear in 
pairs. As the hopping increases, the distance between the nodes in a pair decreases. When 
(H) cannot be satisfied the nodes are aligned. The situation is different for 0o = 7r /2 where a 
continuous line of nodes is obtained. It can be seen that as the hopping increases more values 
of satisfy (0) and when t > A there are nodes for any 0. This behavior is shown in Figs. 
2c and 2d. The continuous line of nodes implies the existence of a finite density of states at 
E = 0. This is only observed for values of 0o very close to tt/2. In a superconductor with 
tetragonal symmetry values of O > tt /4 correspond to non-equilibrium situations. However, 
they can appear if a current is applied, as well as in orthorhombic superconductors with non- 
equivalent axes, where similar features are expected. 

The density of states shows a characteristic multi-peak structure due to van Hove singu- 
larities, which are associated to saddle points of E^ given by eq (0). It is worthy to note that 
we have started from a which presents no saddle points. As can be seen in Fig. 3, the 
width of the peaks increases with the hopping reflecting the splitting of the Fermi surface. 
For high energies the density of states is a constant as in normal metals in 2D. 

In general there are localized states at the grain boundary. Notice that we use localized 
states in the sense that for a given momentum there exist excitations with an energy Ek 
lower than the corresponding excitation energy in absence of grain boundary (0 O = 0). 

E kl = ^/(&±t) 2 +| A~Tp (8) 

where 7 labels the index of the plane. The region of the Brillouin zone which presents 
localized states is displayed in Fig. 4. These states are mainly found inside the splitted Fermi 
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surface (£& = ±t) for the unfrustrated case and close to the values of which correspond 
to the nodes. The area of this region increases with t and 4> . This kind of grain boundary 
in a 3D sample could be described as two coupled semiinfmite media. The localized states 
could not propagate into the bulk, but they would only exist at the boundary. The effect 
of the semiinfmite media can be introduced in terms of a self-energy. As a consequence, the 
extended states will have a finite mean lifetime. Note that, in a real sample, the planes are 
coupled in groups. In some cases it is only necessary to consider the coupling inside the 
groups. Thus, our results could be directly applied to double-layered systems. 

In conclusion, we have studied a new type of grain boundary, which can appear in gran- 
ular superconductors and, in general, during the growth process. The excitation spectrum 
displays new interesting features which seem to be common to other grain boundaries. Our 
model system can be found in any superconductor with d-wave symmetry and not only in 



the ones with orthorhombic structures as the widely studied twin boundaries [p^|15| , piH 2^] 
in YBCO. The simplicity of our model allows us to explore the whole spectrum. In more 
complex systems where not only two planes but two semiinfinite media are coupled, the spec- 
trum is expected to display similar properties at the boundary. The quasiparticle spectrum 
has nodes such that the positions of the latter depend on the coupling constant between 
the planes and on the relative angle of the order parameters. The multi-peak structure 
of the density of states reflects the saddle points of the quasiparticle spectrum, while the 
enhancement at low energies is due to the existence of localized states. 
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FIGURES 

FIG. 1. Sketch of our model system 

FIG. 2. The nodes of the spectrum for O = f , t = 0.4A in Fig. 2a and t = 0.8 A in Fig. 26, 
and for </>o = § , t = 0.8A in Fig. 2c and t = 1.2A in Fig. 2d, are plotted. (j>o = § displays an 
unusual behavior with a continuum line of nodes in contrast to eight nodes found in almost all the 
other cases. The influence of the hopping can also be seen. For <£ = f the nodes in a pair tend 
to be closer as the hopping increases (see text). For </>o = § as the hopping grows an increasing 
number of nodes is obtained. Here, k is given in units of a -1 , where a is the lattice constant. 
Finally, = 40A and n = 10A. 

FIG. 3. The density of states for 0o = f , t = 0.8A and t = 2 A is plotted. The enhancement 
of the density of states at low energies increases with the hopping. The hopping also influences the 
width and number of peaks which appear. 

FIG. 4. The figure shows the region of the Brillouin zone (BZ) which displays localized states 
for the case <fo = f an d t = 0.8A. Units of k and values at the boundaries of BZ are the same as 
those given in the caption of Fig. 2 
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Fig. 1 
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